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C^ . Abstract. - Westudy the current-phase relation of the superconductor/normal/superconductor 

(SNS) junction imbedded in a superconducting loop. Considering the current conservation and 
free energy minimum conditions, we obtain the persistent currents of the SNS loop. At finite 
' O ' temperature we can explain the experimentally observed highly non-sinusoidal currents which 

C ' have maxima near the zero external flux. 

o 

. !^ , Introduction. - Physically interesting superconducting loops with junctions are the su- 

^^ ' perconducting quantum interference device (SQUID) and the superconductor/normal/superconductor 
Tij" , loop with long normal sector (SNiongS loop). Both are mesoscopic size and yield persistent 
^«0 ' currents when a magnetic flux threads the loop. The thickness of the Josephson junctions in 
~^ , SQUID loop is much smaller than the superconducting coherence length, but the thickness 
pi ' of normal segment d of SNiongS loop can be larger than that. The persistent current in the 
former system flows by tunneling the Josephson junction, while that in the latter by the long 
C^ ' range proximity effect. Since these two loops are very different in nature, we can anticipate 
a different current-phase relation for SNiongS loop, but it was not explicitly demonstrated in 
previous studies for the superconductor/normal/superconductor hybrid junction [1-3] and for 
' ^ I the SNiongS loop [4,5] in connection with the Andreev reflection process [6]. 

The current described by the motion of a pair of electron and hole in the normal sector 
^ , changes into that by a Cooper pair in the superconducting sector. There is an intermediate 
region near the edge of superconductor, where the current is described by quasiparticles and 
quasiholes.Then the current carried by the electrons and the holes in the normal sector should 
be the same as that carried by quasiparticles and quasiholes in the intermediate region. This 
current conservation condition can be satisfied by considering that two wave vectors of the 
quasiparticle and the quasihole can be different from each other like those of electron and hole 
in the normal sector. In solving the Bogoliubov-de Gennes (BdG) equation, we consider the 
average energy of a pair of particles is a dynamic variable rather than a constant chemical 
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potential. Thus the values of the dynamic variables are determined by minimizing the free 
energy of the SNiongS loop. 

Recently an experiment for the high- Tc supercondutor (HTS) junction with interlayer much 
thicker than the superconducting correlation length has been reported [7], where the HTS 
junction was incorporated in a superconducting loop with threading external magnetic flux 
$ext. In this single-junction interferometer experiment, the current-phase relation shows the 
highly non-sinusoidal behaviors so that the slope of the current for zero external flux becomes 
larger than that for (f>oxt — $0/2 with the superconducting unit flux quantum $0 = /i/2e as 
temperature goes down. The tunneling currents are supposed to take place through the long 
range proximity effect across the thick PrBa2Cu307 (PBCO) interlayer. Since the CuO chains 
in the PBCO layer are metallic at low temperature, the above junction can be considered as 
a SNiongS junction incorporated in a superconduction loop. In this study we calculate the 
persistent current of the above SNiongS loop considering the current conservation and free 
energy minimum conditions and show that the experimentally observed non-sinusoidal type 
current can emerge through the long-range proximity effect in the SNiongS loop. 

Andreev reflections in a superconductor/normal/superconductor loop. - The quasiparti- 
cles in an SNS loop with threading external magnetic flux can be described by the Bogoliubov- 
de Gennes (BdG) equation 

H^-i A{z) \ f u{z)\ _ f uiz)\ 

where Hq = {—ihd/dz — eA/cf' /2me with the electron mass me, A = ^cxt/L is the vector 
potential with the circumference of loop L, and A(z) is the pair potential as a function of 
spatial coordinate z. 

For a superconducting loop without junction threaded by an Aharonov-Bohm flux, the 
quasiparticle wave functions u{z) and v{z) contain extra factor due to gauge invariance such as 
u{z) = Ce'(A=+'^//i)^ and v{z) = Ce^(^h-7r//L)2:^ ^^^^.^ j = $ext/$o- Uniform flow of persis- 
tent current is derived by the BdG equation with the pair potential, A(z) = /\g^{^i~^h+-^f /l)z 

If the superconducting loop is interrupted by a normal sector as shown in Fig. ^ the wave 
function ^ = {u{z) v{z))'^ of a pair of electron and hole in normal sector {0 < z < d) and 
of a pair of quasiparticle and quasiholc in the intermediate region of superconducting sector 
(d < z and z < L) is given by 



u(z) 
v{z) 



(^Cgz(A+.-j/).-^^j id<z) (2) 

JJ^iiX-H-ff)z+i,i J {z<L), 



where r] is the phase shift due to Andreev reflection at the NS interface [4,5,8] and, when a 
pair of quasiparticles passes through the NS interface, each quasiparticle acquires the addi- 
tional phase equal to rj. Since we do not make any assumption about the sign of the wave 
vectors fco and fci, this wave function can describe both the excitations moving clockwise and 
counterclockwise. Here a notable point is that the wave vectors of the quasiparticle and the 
quasihole in the previous work [2, 4, 5] are set to be the same. It is, however, natural to 
discriminate the wave vectors like those of particles /cq and fci in the normal sector in order to 
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Fig. 1 - A superconductor/normal/superconductor loop with length L. k and A denote the wave 
vector of a electron in normal sector and a quasiparticle in superconducting sector, respectively. 

satisfy the current conservation condition. Therefore, we introduce different wave vectors, Ag 
and A;i, for the quasiparticle and the quasihole in the intermediate region of superconducting 
sector. 

This equation can be solved easily for the normal sector in which A{z) = 0, and yields the 
relations for ^ and E such that 



E-- 



Anie 
Am, 



-(/Cq — ki). 



(3) 
(4) 



In the intermediate region of superconducting sector, however, the BdG equation must be 
solved with the pair potential 



^^' ^ 1 ^gi(A_,-A_;.+27r//L)2 (^^ <L). 



(5) 



Then the BdG equation for z > d becomes h^\\^/2m(, — ^ + Ae^^"' = E and Ae^"' — 
h^\\i-^/2rae + ^ — E. Representing A+e and A+h such as A+e = ^o + i^'o and A+^i = Ai + iX[ 
we get an expression for ^ and E, 



e 



E = 



4m, 
Am, 



iXl + Xl)il-a) 



(Ag - Xl){l + a) + Acos277 



(6) 
(7) 



with a = {meAsm2ri/h^XoXi)^, Xq = We A sin 2?7/fi,^Ao and X[ — m,eAsm2ri/h^Xi. For 

z < L, we can also solve the BdG equation and find that A_e = Aq — iXq and A_;j = Ai — iX[. 

The phase matching conditions for the wave function of Eq. ^ at z — d and z = L are 

given bv Ae^*'"'^"*'^'''*"*'' = J^Q^f^id-iXid+ij-i _ ^ g^^^([^Q^ii.^o+Trf/L)L+iri _ ^g-j(Ai-7r//L)L-ir; _ 

D, where A, B, C and D represent the phase part of the coefficients. A, B, C and D, respec- 
tively. The condition for the existence of a solution leads to the boundary condition 



ko - ki 



¥)- 



Ao - Ai + -y!- \ {L-d)-AT] = 27m. 



(8) 
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When an electron becomes Andreev-refiected at the superconductor/nomal interface, the 
transmitted quasiparticle pair obtains the additional phase 277. The phase Ai] in the boundary 
condition in Eq. ^ is the sum of two phase changes 2?/ due to the Andreev reflections at 
each interface. 

The condition of current conservation at NS interface will be given by using the represen- 
tation of the flux j = —{l/2m)[(j){—ihd/dz — eA/c)(j>* — (j)*{—ihd/dz — eA/c)(t>\. In the wave 
function '5 — {u{z) v{z)Y' in Eq. 0, u{z) is the electronlike wave function and v{z) is the 
holelike wave function which is the complex conjugate of the electronlike wave function. The 
current representation for a pair of electrons thus should be obtained with the wave function 
4> = {u{z) v*{z)Y' and we can get the relation 



fci = Aq — Ai = 2g. 



(9) 



The Cooper pairs in superconducting sector then should carry the current / = nc'2.eh2q/mc 
with the density of Cooper pairs Uc and ttic — 2me- 

Persistent currents of a superconductor/normal/superconductor loop. - When we calcu- 
late the current through the SNiongS junction, we may be able to set the average energy ^ 
of the particle and the hole in Eq. Q equal to the chemical potential. But, in the SNiongS 
loop, S, in Eq. (^J need not be the constant chemical potential. For example, we consider 
a simple normal loop with a threading flux. The average energy of two particles S, at Fermi 
level is different from the chemical potential /i such that ^ = {l/2){h^ /2me){k'^ + k^) = 
(l/2)(;iV2me)((fcF - 7r//L)2 -|- [kp + T^I/Lf) + (h^ I2m^)k\ = /i. Since fco and fci depend 
on the dynamic variable r\ as well as / in the SNiongS loop, ^ cannot be set as a constant 
chemical potential but should also be a dynamic variable to be determined. 

Since an extra variable ^ is introduced, we need one more independent relation. That 
is given by the requirement of free energy minimum [9]. Since the intermediate region of 
superconducting sector is so thin that we neglect the energy of this region in the free energy 
expression and consider only the energy of Cooper pairs. In superconducting sector the Cooper 
pairs carry the persistent current corresponding to the Cooper pair wave vector 2q = \q — Ai 
and the energy of a Cooper pair can be written as {h^ /2mc){2q)'^ with rric — 2m^. Therefore 
the total free energy per particle f/tot can be written as 



2C/tot = 



2?Tip 



(fco^ + fcf)- + ^^(Ao-Ai)^(l- 



(10) 



Using Eqs. Q and ^, the total free energy can be represented in Aq and Ai. 
From the condition dUtot/drj — 0, we obtained the equation. 



(X +\) + (\ A )^ I (2"^eAsin2r?)^ l 
(Ao + Ai) + (Ao - Ai)^ + ^^p 

(Ao + Ai) - (Ao - Ai)- + ^^, 



dXp 
drj 
dXi 
drj 



Xl + Xf f2meA 



2\2 



2A0A1 



ft2 



sin 477. 



(11) 



Here dXo/dr/ and dXi/dr/ can be obtained by differentiating Eqs. Q, (©, ©, 0, I® and © 
such as dAo/d77 = [jf'+{8Ai(afco + fci)+4Aorsin27?}/L]/2?anddAi/d?7 = [T+{8XQ{ako + ki)- 
4Airsin277}/i]/r', where J^ = 2AoAircos2r/-2A(fco-Kfci) sin277-8fcofci/£ and P ee -2(Ao- 
Xi){ako + fci) + (Ao + Ai)rsin2?7 with F = (fcoAg - fciAf)A2 sin2?7/2AgA? and A = 2meA/h'^. 
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Fig. 2 - (a) Excitation spectrum of the SNiongS loop. The sohd (dashed) hne for (^3 < ((/? > 0) shows 
the ground state energy, where y = 27r/. (b) Persistent currents corresponding to the energy levels 

in (a). Here we set Lv Ao = 10 and A/Ao = 1. 



Solving the coupled equations ©, Q, ©, Q, ©, ® and l(TT)) numerically, we obtain 
the energy levels corresponding to the two solutions of the BdG equation in Eq. (|T|l as a 
function of external flux ip = 2tiJ — 27r$ext/'&o as shown in Fig. Efa), where the solid and 
the dashed lines denote the lower levels and the other lines the higher levels. The ground 
state corresponds to the solid (dashed) line for (^ < (v? > 0). The persistent currents 
in Fig. |2Ib) are represented by the same lines as those of the corresponding states in Fig. 

121a). The persistent current, / = nc2eh2q/rnc, can be written by / = /o2g/\/ Aq, where 
Jo = ric^ehy Kq / rric and Aq = 2r7ZeAo//i^ is a gap potential chosen arbitrary. 

Since the Cooper pairs in the superconductor are in the coherent condensate state, each 
Cooper pair carries the same superconducting current. A Cooper pair changes into a pair of 
normal electrons in the normal sector via a pair of quasiparticles in the intermediate region. 
Since the current in the loop should be conserved, the current equal to the macroscopic 
persistent current in the superconducting sector flows in the normal sector. In Fig. El^b), we 
can see several persistent currents, of which the persistent current at ground state corresponds 
to the solid line ior ip <Q and the dashed line for (^a > 0, the saw-tooth type current. 

Recently a single-junction interferometer experiment was done on HTS junction [7] (YBa2Cu307_x 
/PrBa2Cu307/YBa2Cu307_x) which is incorporated into a superconducting loop with pene- 
trating magnetic flux. Since there is no misorientation angle between the two d-wave super- 
conductors across the inter layer in this experiment, the phase difference across the interlayer 
can be brought about only by the threading magnetic flux. They obtained the current-phase 
relation at finite temperature, where the slope for zero external flux is larger than that for 
the external flux ip ~ tt. 

In this experiment, the thickness of the PBCO interlayer of the junction is as long as 
hundreds A which is much larger than the high- Tc superconducting correlation length by an 
order of magnitude and thus the Cooper pairs cannot directly tunnel the junction. These long 
range proximity effects have been observed in many experiments on HTS junctions [10-13] 
and considered as a characteristic of SNjongS junctions. Thus the experimental results may 
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Fig. 3 - Current-phase relation at finite temperatures. As the temperature goes down, the maximum 
approaches the point </? = 



be explained by calculating the currents of the SNiongS loop as a function of the threading 
external flux. 

Here we can raise a question that whether it is a real proximity effect, i.e., there are 
filaments of pinholes or microshots so that the Cooper pairs can be transported by the resonant 
tunneling through localized states in filaments. Recently an experiment on the trilayer HTS 
junction [11] has been reported, where they synthesized atomically smooth films of HTS and 
uniform trilayer junction so that they conclude that the long-range proximity effect do not 
originate from the resonant tunneling through the energy-aligned states in the microshots, 
but it is an intrinsic property of the interlayer. In fact, the supercurrents are known to flow 
via metallic CuO chains in the PBCO interlayer [14,15]. 

From the currents in Fig. [2Ib) corresponding to the states in Fig. [2Ia) we can obtain the 
persistent currents of the SNjongS loop at finite temperatures. At finite temperature, since 
the current state with energy e has the probability proportional to e"'^/''^^, we can obtain 
the thermally averaged persistent currents from those current states in Fig. [21 as shown in 
Fig. 121 We can observe highly non-sinusoidal currents in Fig. O where the amplitude maxima 
approach the point 93 = and the slope of current at (/? = becomes larger than that at (/3 = tt 
as temperature goes down. The current takes the sinusoidal form only after the temperature 
goes up such that fc^T/A « 0.4 which corresponds to the temperature near Tc- 

Actually the similar behavior can also be seen in the loop interrupted by a grain boundary 
Josephson junction not by an SNiongS junction, which however comes from the phase difference 
across the grain boundary Josephson junction due to the misorientation angle between two 
d-wave superconductors [16]. In SNiongS junctions the wave functions acquire the Andreev re- 
flection phase shift r\ which appears in the boundary condition of Eq. . The Andreev phase 
shift takes the place of the misorientation angle of the grain boundary Josephson junction 
and thus results in the non-sinusoidal current-phase relation as shown in the manuscript. In 
Ref. [7] the authors explained their own experimental results by assuming thermal fluctuations 
which induced the highly non-sinusoidal current-phase relation at low temperatures. At higher 
temperatures near the critical temperature they can recover the sinusoidal current-phase re- 
lations. In the present study, however, we can show that the non-sinusoidal current-phase 
relation emerges naturally at low temperature without assuming thermal fluctuations. 
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Summary. - We have studied the current-phase relation of the SNiongS loop with thread- 
ing magnetic flux. The Cooper pairs in the superconducting sector become the electron pairs 
in the normal sector via the quasiparticle states in the intermediate region. The net current in 
the SNiongS loop should be preserved and thus we introduce different wave vectors of the pair 
of quasiparticles in order to satisfy the current conservation condition. Furthermore, since 
the average energy of a pair of particles should be a dynamic variable, we also consider the 
free energy minimum condition. We obtained the persistent currents of the ground and ex- 
cited states and found that we can explain the experimentally observed highly non-sinusoidal 
current-phase relation with maxima near the zero external flux at low temperature. 
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